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Abstract 

We study bifurcations of homoclinic orbits to hyperbolic saddle 
equilibria in a class of four-dimensional systems which may be Hamil- 
tonian or not. Only one parameter is enough to treat these types of 
bifurcations in Hamiltonian systems but two parameters are needed 
in general systems. We apply a version of Melnikov's method due to 
Gruendler to obtain saddle-node and pitchfork types of bifurcation 
results for homoclinic orbits. Furthermore we prove that if these bi- 
furcations occur, then the variational equations around the homoclinic 
orbits are integrable in the meaning of differential Galois theory un- 
der the assumption that the homoclinic orbits lie on analytic invariant 
manifolds. We illustrate our theories with an example which arises as 
stationary states of coupled real Ginzburg-Landau partial differential 
equations, and demonstrate the theoretical results by numerical ones. 
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1 Introduction 

Bifurcations of homoclinic orbits in ordinary differential equations have been 
studied in numerous articles over the past decades. They also arise as bi- 
furcations of solitons or pulses in partial differential equations (PDEs), and 
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Figure 1: Assumption (A2) 



have attracted much attention even in the fields of PDEs and nonlinear waves 
(see, e.g., Section 2 of [28]). 

Many of researches on this topic are related to bifurcations from single- 
bump homoclinic orbits to multi-bump ones. See Section 5.2.4 of [28] for 
a concise and useful review of the references. As an exceptional study, 
Knobloch [15] showed that a saddle-node type of bifurcation for homoclinic 
orbits to hyperbolic saddles can occur in reversible and conservative systems 
with one parameter when the homoclinic orbits are non-transversal, i.e., the 
tangent spaces to the stable and unstable manifolds of the saddles along 
them have two-dimensional intersections, using so-called Lin's method [18]. 
We can also use a version of Melnikov's method [HI [19] due to Gruendler [TO] , 
who studied the persistence of homoclinic orbits, to treat such bifurcations 
in general higher-dimensional systems with two parameters. We note that 
only one parameter is enough to treat these bifurcations for reversible and 
conservative systems but two parameters are needed in general systems. A 
pitchfork type of bifurcation for homoclinic orbits to saddle-centers was also 
detected in reversible systems [31], using an idea similar to that of Melnikov's 
method. 

Here we are interested in the latter, saddle-node and pitchfork types of 
bifurcations for homoclinic orbits to hyperbolic saddles in systems of the 
form 

x = f(x;fi), iGl", /i6l m , (1) 

where / : M. n x W m — > M. n is analytic, /i is a parameter (vector), n is a positive 
integer and m = 1 or 2. Note that Eq. (pQ) may be Hamiltonian or not. We 



2 



consider the case of n = 4 and make the following assumptions. 

(Al) The origin x = is a hyperbolic saddle equilibrium in ([T]) at // = 0, 
such that D x /(0; 0) has four real eigenvalues, Ai < A 2 < < A 3 < A 4 . 

(A2) At /jL — the hyperbolic saddle x = has a homoclinic orbit x h (t). 
Moreover, there exists a two-dimensional analytic invariant manifold 
^# containing x = and x h (t) (see Fig. [T]). 

From (Al) we can also assume that the origin x = is still a hyperbolic saddle 
near /i = under some change of coordinates. Assumption (A2) implies that 
the complexification of the vector field ([I]) on ^# with /x = has a complex 
separatrix T whose real part is the real homoclinic orbit x h (t). In particular, 
the existence of the invariant manifold ^ is essential in Theorem 11.11 below. 
We also remark that as well known [T2 ] I2T ] 129] . there exists a center manifold 
for x h (t) under a more general condition, but it is not analytic in general. 

As seen in [TU] US ESj, variational equations (VEs) play an important 
role in persistence and bifurcation problems of homoclinic orbits. The VE of 
(II]) around the homoclinic orbit at /i = is given by 

e = D x /(x h (t);0)£, (2) 

where £ G M n . We easily see that £ = x h (t) is a bounded solution of (T5]) 
tending to zero exponentially as t — > ±oo. We consider the following condi- 
tion. 

1. The VE 02]) has another independent bounded solution. 

In Section 2, using Gruendler's version of Melnikov's method [TU] instead of 
Lin's method [18], we prove the following results for ([1]) with n > 4: 

(i) If condition (C) holds, then a saddle-node or pitchfork bifurcation of 
homoclinic orbits occurs under some nondegenerate condition; 

(ii) If condition (C) does not hold, then these bifurcations do not occur. 

Thus, condition (C) provides a criterion for saddle-node and pitchfork bifur- 
cations of homoclinic orbits in ([1]). 

On the other hand, the problem of integrability is especially important in 
dynamical systems, along with bifurcations and chaos. In particular, it has 
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been extensively studied for two classes of dynamical systems, linear differ- 
ential equations and Hamiltonian systems. See, e.g., [21 [30] for the notions 
of integrability in these two classes. The former is also briefly described in 
Section 3. Differential Galois theory is one for the former class and provides 
algorithms to solve integrable linear differential equations. Morales- Ruiz and 
Ramis [221 [21] applied the differential Galois theory and uncovered a relation 
between integrability of both classes: If a Hamiltonian system is integrable 
near its particular solution in the context of complex Hamiltonian systems 
(i.e., in the sense of Liouville with meromorphic first integrals), then the re- 
lated VE is integrable in the context of linear differential equations (i.e., in 
the meaning of differential Galois theory). Their result was further extended 
to higher-order VEs [2S] . An implication of their result for the occurrence of 
chaotic dynamics in two-degree-of-freedom Hamiltonian systems with saddle- 
centers was also discussed in [231 [33] . 

In this series of articles, we give another application of the differential 
Galois theory and clarify a connection of saddle-node and pitchfork bifurca- 
tions of homoclinic orbits in ([T]) with the integrability of the VE 02]) in the 
meaning of differential Galois theory. More precisely, we prove the following 
theorem (The proof is given in Section 4). 

Theorem 1.1. Suppose that condition (C) holds. Then the VE §2§ has a tri- 
angularizable differential Galois group, when regarded as a complex differen- 
tial equation with meromorphic coefficients in a desingularized neighborhood 
Tioc of the homoclinic orbit x h (t) in C 4 . 

This theorem means that under condition (C) the VE (T5]) is integrable in 
the meaning of differential Galois theory. See Sections 3 and 4 for more details 
on the statement. Such a relation between integrability in the differential 
Galois setting and the existence of bounded solutions for linear differential 
equations was empirically noticed in the context of exactly solvable potentials 
in one dimensional Schrodinger equations by Acosta-Humanez et al. [I]. Here 
we concentrate ourselves on the case of general or Hamiltonian systems with 
hyperbolic saddle equilibria. The case of reversible systems or saddle-center 
equilibria will be treated in part II. 

The outline of this paper is as follows: In Section 2 we extend Gru- 
endler's version of Melnikov's method [10] and present bifurcation results for 
homoclinic orbits not only in Hamiltonian systems but also in more general 
systems. These results are easily obtained from the result of [10] but many of 
them are not found in other references, to the authors' knowledge. We give 
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necessary information on differential Galois theory in Section 3 and state a 
proof of Theorem 11.11 in Section 4. In the proof, under assumption (A2), 
we reduce integrability of the VE (J2]) to that of a two-dimensional linear 
system which consists of components normal to the invariant manifold 
and is called the normal variational equation (NVE). Finally, in Section 5, 
we illustrate our theories with an example which arises as stationary states 
of coupled real Ginzburg-Landau PDEs. The theoretical results are also 
demonstrated by numerical ones using a numerical computation tool called 
AUT097 [7]. 

2 Melnikov analyses 

In this section, we extend Gruendler's version of Melnikov's method [10] and 
discuss the persistence and bifurcations of homoclinic orbits in ([TJ). 

2.1 Extension of Melnikov's method 

We consider a more general situation of n > 2 in ([TJ and assume the follow- 
ings: 

(Ml) The origin x = is a hyperbolic saddle in flTJ at // = such that 
Da;/(0; 0) has n s and n u eigenvalues with negative and positive real 
parts, respectively, where n s + n u = n. 

(M2) The equilibrium x = has a homoclinic orbit x h (t) at fi = 0. 

(M3) The VE (J2J) has just n independent bounded solutions, £ = (pi(t) = 
x h (t),ip 2 (t), . . • ,<£>„„(*). 

Assumptions (Ml) and (M2) mean that the stable and unstable manifolds 
of the origin at ji = 0, Wq(0) and Wq^O), are of dimension n s and n u , 
respectively, and follow from (Al) and (A2) for n = 4 with n s = n u = 2. By 
assuming (M3), we have n s ,n u > n and 

dhn(T x WZ(0) n T x W£(0)) = n 

along the homoclinic orbit x h (t). Note that the existence of an analytic 
invariant manifold containing the equilibrium x = and homoclinic orbit 
x h (t) is not required. 

Using Theorem 1 of [TU], we obtain the following lemma immediately. 
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Lemma 2.1. There exists a fundamental matrix = (<pi(t), . . . , <p n {t)) of 
(TJP snc/i that for nonnegative integers kj, j = 1, . . . , n, and a permutation a 
on n symbols {1, . . . , n}, 

<p j (t)t~ k ie- ReCx i )t = 0(1) as t -> +oo, 
^■(tJt-^OJe-^W))* = as t -> -oo, 

where Xj, j = 1, . . . , n, denote the eigenvalues o/D x /(0; 0) repeated according 
to algebraic multiplicity, and Re(Aj) is negative if j < n s and positive if 
J > n u . 

By assumption (M3) and Lemma 12.11 we have 



lim if At) = for j = 1, . . . , n ; 

t— >±oo 

lim ^ jit) = 0, lim <Pj(t) = oo for j = no + 1, . . . , n s 

t— >+oo t— >— oo 



(3) 



and assume that 



lim <Pj(t) = oo for j = n s + 1, . . . , n s + n ; 

t— >±oo 

lim <fj(t) = oo, lim <fj(t) =0 for j — n s + n + 1, . . . , n. 

t— >oc t— oo 



(4) 



Define ipj{t) for each j = 1, . . . , n by 

{ip jit), ip k (t)) =6 jk , 

where (5^ is Kronecker's delta, and (£, rj) represents the inner product of 
£,7/ G M n . It immediately follows from fl3]) and fl4]) that 

lim V'jW = oo for j = 1, . . . , n ; 

t— >±oo 

lim ^j(i) = oo, lim ipj(t) = for j = no + 1, . . . , n s ; 

t— > + 00 t— > — OO /r\ 

lim ipj(t) = for j = n s + 1, . . . , n s + no] 

t— >±oo 

lim ^>,-(i) = 0, lim ipAt) = oo for j = n s + no + 1, . . . , n. 

t-^oo t— >— oo 

The functions ipj(t), j = 1, . . . ,n, can be obtained by the formula *$?(t) = 
where *ff(t) = (ip 1 (t),. . . , ip n (t)) and $*(t) is the transpose matrix 
of $(£). Moreover, ^(t) is a fundamental matrix to the adjoint equation 

Z = -D x f(x h (t);0)*Z (6) 
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since ty*(t)$(t) = id n with id n the nxn identity matrix, so that W*(t)$(t) + 
y*(t)&(t) = 0, i.e., 

if(t) = -(**(t)$(t)$-\t)y = -D x f(x h (ty,oyv(t). 

Now we look for a homoclinic orbit of the form 

no— 1 

x = x h (t) + a m+i(t) + ^(VW+W) (7) 

3=1 

in ([[]) with (j ^ 0, where a = (oti, . . . , a no _i). Here the ^(a)-terms are 
ignored in (J7J) if n = 1. 

Let k be a positive real number such that 

k< J |Re(Aj)|, j = l,...,n, 

and define two Banach spaces 

iT = {z g C°(1R, lR n ) | sup |,z(t)|e K| * 1 < oo}, 

t 

= { z 6 C^R, R n ) | sup Iz^le" 1 '*' < oo,swp\z(t)\e K ^ < oo}, 

where the maximum of the suprema is taken as a norm of each spaces. The 
following lemma is frequently used in the rest of this section. 

Lemma 2.2. Let z G Sf 1 . Then we have 

(i> ne+j {t), z{t) - B x f(x h (t); 0)z(t)) dt = 0, j = 1, . . . , no. 

Proof. The proof is found dividedly in [TO]. Here it is given for the reader's 
convenience. 

We easily see that for z G iF° 

|(VWt),z(*)>| <K \\z\\e~ Kt , j = l,...,n , 

by Lemma T2. II and (jSJ), where i^o is some constant independent of z. On the 
other hand, 

(i> ns+j (t),z(t) -D x f(x h (t);0)z(t)) 

= (i> ns+j (t),z(t)) - (D x f(x h (t);0)*i> n3+j (t),z(t)) 

= {^ ns+j {t),z{t)) + (iP ns+J (t),z(t)) = A(^ s+J (t),z(t)> 

since ip ns+ j{t) is a solution (jSj). Thus, we obtain the result. □ 
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From Theorem 4 of [TO] we also have the following lemma. 
Lemma 2.3. The nonhomogeneous VE, 

£ = D x f(x h (t);0)Z + V (t) (8) 
with n G iF° , has a solution in iF 1 if and only if 

(,/> nB+j (t),7i{t))dt = 0, J = l,-..,n . (9) 

Moreover, if condition (TJ|) holds, then there exists a unique solution to (TJj] 
satisfying (^j(O), f (0)) = /or j = 1, . . . , n . 

Let 

^ = {ze^ l \ (^(0), z(0)> = 0, j = 1, . . . , n Q } C & 1 

As the kernel of a continuous linear map 3f l is closed and hence it is also a 
Banach space. We define a differentiable function F : 3f 

as 

no-l 

F(z; a,/i) =x h (t) + + a j L Pj+i{ t ) 

3=1 

/ "o-l \ 

- / (* h (*) + + £ -"J • ( 10 ) 

A solution z G iF 1 to 

F(z;a,/i)=0 (11) 

for (a,n) fixed gives a homoclinic orbit to the origin. 
Let q : R — >■ M be a continuous function satisfying 

/oo 
g(t)dt = 1. 
•00 

Define a projection IT : — >• iF° by 

n^(t) = g(t) ^ (j (^n s +i(r), z(r)) drj ^ S+J (t). 



Condition (II ip is divided into two parts: 

(id - U)F(z;a,ii) = 



(12) 



and 



[[^(z; a, fi) = 0, 
where "id" represents the identity. Obviously, 

(id - n)F(0; 0, 0) = 0. 

We easily see that 



(13) 



(14) 



/ {^ ns+j (t),(id-U)z(t))dt = 0, j = l, 



...,n , 



(15) 



j —00 



for z G Hence, by Lemma I2.3[ Eq. ([8]) has a unique solution in 

if t] is involved in the range of (id — II), ^(id — II). This means that for 
r] G ^(id — II) there is a unique function £(t) G S"\ such that 



Using this fact and f[T4]) . we apply the implicit function theorem (e.g., The- 
orem 2.3 of [6]) to show that there are a neighborhood U of (a,/x) = (0,0) 
and a differentiable function z : U —> iF 1 such that z(0, 0) = and 



Fj(a,fi)= (ip ns+j {t),F(z{a,fi){t);a,fi))dt, j = l,...,n . (17) 



If F(a, fx) = (Fi(a, //),..., F no (a, /x)) = 0, then z = z(a, /x) satisfies ( IT^|) and 
f ll3p . Thus, we can prove the following theorem as in Theorem 5 of [TU]. 

Theorem 2.4. In a neighborhood U of (ct,/x) = (0,0) , Eq. ([TJ /ias a unique 
homoclinic orbit of the form (J7J) to i/ie origin x = if F(a, /x) = 0. 

Henceforth we apply Theorem 12.41 to prove persistence and bifurcation 
theorems for homoclinic orbits in §Q) with n > 2. The first two results 
(Theorems 12.51 and I2.7P are essentially the the same as Corollaries 7 and 8 
of [10] but they are included with proofs for our unified presentation and 
the reader's convenience. We will also need information on higher-order 
derivatives of F up to third-order at (a, /x) = (0, 0), some of which were not 
required in [TO] . 



D 2 (id-n)F(0;0,0)£ = ?7. 



(id — H)F(z(a, /x); a, /x) = 



(16) 



for (a, /x) G U. 
Let 
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Figure 2: Persistence of homoclinic orbits 



2.2 General case 

Let m = 2 and assume that no = 1, which means that condition (C) does 
not hold when n — 4. Then F is independent of a. Define two constants 
an, ai2 € K as 

/oo 
(^ ns+1 (t),D Mj /(x h (t);0))dt. (18) 
-oo 

Theorem 2.5. Under assumptions (M1)-(M3) let uq = 1 and suppose that 
Oi = (011,012) 7^ (0,0). T/ien /or some open interval I including 0, there 
exists a differentiable function <f)j : I — > R wii/i 0^(0) = 0, j = 1 or swc/i 
i/iai a homoclinic orbit exists for fix = 4>\{iAi) or fi 2 = ^2(^1) r see [IP- 

Proof. We simply write 2 = z(fi). For /i 6 I 2 we compute (fl71) as 
AG") 

/oo 
- D x /(x h (t); 0)z(t) - BJ(x\t); 0)fj,)dt + €?(|/i| 2 ) 
-00 

= -ai/i+ ^(|/i| 2 ). 

Here we used the fact that z(0, 0) = and Lemma 12.21 Applying Theo- 
rem [231 we obtain the result. □ 

Remark 2.6. Theorem \2.5\ implies that if condition (C) does not hold for 
n = 4, then the homoclinic orbit x h (t) persists, i.e., no bifurcation occurs, 
since no = 1, as stated in Section 1. 
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Figure 3: Saddle- node bifurcation of homoclinic orbits 



We next assume that n = 2, which means that condition (C) holds when 
n — 4. For /i G M 2 define constant vectors ctj G M 2 and constants bj G M, 
j = 1,2, as 



Note that ai = (an,a 12 ), where ay, j = 1,2, are given by (|T8l) . as in the 
statement of Theorem 12. 51 

Theorem 2.7. Under assumptions (M1)-(M3), let n$ = 2 and suppose that 
the 2x2 matrix (a^a^) is nonsingular and (61,62) 7^ (0,0). Then for some 
open interval I including 0, there exists a differentiable function : / — > M 2 
with 0(0) = 0, 0'(O) = and 0"(O) 7^ 0, such that a homoclinic orbit exists 
for fi = 0(a), i.e., a saddle-node bifurcation of homoclinic orbits occurs at 
H = (see Fig. 

Proof. Differentiating (I16p with respect to a and using Lemma 12. 2\ we have 



at (a, 11) = (0,0), i.e., D a z(0;0)(t) is a solution of ©, so that D a f(0;0)(t) = 




(19) 



D a (id - U)F(z; 0; 0) = -D Q z - B x f(x\t); Q)D a z = 



11 



by Lemma 12.31 Using this fact and Lemma 12.2} we compute (I17p as 

/oo 
(^ 8+i (t),-D M /(x h (t);0> 
-oo 

- \a 2 V 2 J{x\t)- 0)(^(*), Mt)))dt + ^(v^+H 1 ) 
= - CLjfi - bjct 2 + ff(^Ja 6 + |/i| 4 ), 

as in the proof of Theorem 12.51 Since F(0, 0) = and |D M F(0,0)| 7^ 0, we 
apply the implicit function theorem to show that there exist an open interval 

1 3 and a different iable function : I — > M 2 such that F((p(a), a) = for 
a € I with 0(0) = 0, 0(0) = and 

This implies the result along with Theorem 12.41 □ 



2.3 Z2-equivalent case 

We next consider the Z^-equivalent case for uq = 2, and assume the following. 

(M4) Eq. ([TJ is ^-equivalent, i.e., there exists an n x n matrix S such that 
S 2 = id n and Sf(x; /i) = f(Sx; /i). 

It follows from assumption (M4) that if x = x(t) is a solution to ([I]), then 
x = Sx(t) is so. We say that the pair x{t) and Sx(t) are S-conjugate if 
7^ Sx(t). See, e.g., Section 7.4 of [17] for more details on Z 2 -equivalent 
systems. In particular, the space R n can be decomposed into a direct sum, 

w n = x + © x~ , 

where S*x = 2 for x G X + and S*x = — x for x G X - . Under assumption (M4), 
if x h (t) G X + for any t G M, then Eq. ([2]) is also Z 2 -equivalent about £ 
since SD x f(x;fi) = D x f(Sx; fi)S in general. Here we need the following 
assumption. 

(M5) For every t G R, x h {t), ?/>„ s+ i(t) G X+ and <^ 2 (t), ^„. s + 2 (t) G X". 

Assumption (M5) also means that ^1 (t) G X + . Moreover, a homoclinic 
orbit of the form ([TJ has an S-conjugate counterpart for a / since it is 
not included in X + . Actually, we cannot apply Theorem 12.71 in this 
follows. 
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Lemma 2.8. Under assumptions (M1)-(M5), we have 

BJ(x h (t); 0)/!, B 2 J(x h (t); 0)(<p 2 (t), G X + 

and 

D,D x f(x\ty, 0)(/i, <p 2 (t)), nlf(x\t)- 0)(<p 2 (t), <p 2 (t), <p 2 (t)) e x- 
for any teM and /i G R 2 . In particular, a 2 = and b 2 = 0. 
Proof. By the Z 2 -equivalence of ([I]) we compute 

D M /(x;0) = SBJ(Sx;0), D 2 x f{x;0)(U) = SB 2 J(Sx; 0)(S£, ££) 

and 

D,D M /(x; 0)(£, /i) = 5 DJV^ 0)(Sf , //), 
D2/(x; 0)(£, £, = SD*/(Sx; 0)(S£, S£, SQ. 

Substituting u = x h (t) and v = (p 2 (t) into the above relations and using the 
condition that x h (t) G X + and ^2^) G X~ , we have 

D A1 /(x h (t);0) = 5D M /(x h (t);0), 

D2/(x h (t);0)(^ 2 (t),^(t)) = SDlf(x h (t);0)( V2 (t),<p 2 (t)) 

and 

D x .D,/(x h (t);0)(^(t),/i) = -5D M /(x h (t);0)(^(t),/i), 

D^(x h (t); 0)fo> 2 (t), ^(f), ^(t)) = -SD 2 J(x\t); 0)& 2 (t), <p 2 {t), <p 2 {t)), 

which implies the statement. □ 

Let £ = and £ Q (t) be, respectively, unique solutions to 

t = D x f(x\t)- 0)£ + (id - n)D M /(x h (t); 0)/i 

and 

£ = B x f(x h (t); 0)£ + ±(id - n)D^(x h (t); 0)(y? 2 (t), <p 2 (t)) 

in iF 1 . It is guaranteed by Lemma 12.31 and ( fl5l) that these solutions exist. 
Denote 

G,(t) = D M /(x h (t);0)/i, 6 (f) = lB 2 J(x h (t);0)(Mt),Mt))- 
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Using the fundamental matrices to (J2]) and $(£) and ^(t), and noting 
that #(t) = we obtain 



f ^(r)(id-n)G M , Q (r)dr + vl>*(0)C Q 



where G M n are constant vectors such that (^(0), £q ) = 0, j = 1,2, 
and £ M ' a (t) are bounded on (—00,00), i.e., one can write £o'° = —$(0)2^'° 
with 5 M ' a = . . . , given by 



/- p—00 



Let 



a 2 /i 



(^•(t), (id - n)G^ Q (t)) dt for j = 3, . . . , n s ; 

{ij)j(t), (id - n)G^ a (i)) dt for j = n s + 3, . . . , n; 

otherwise. 

00 

(^ s+2 (t), |D^(x h (t); 0)(^ 2 (t), ip 2 (t), ip 2 (t)) 

00 

+ D^(x h (t);0)(r(t),^ 2 (t)))dt. 



(20) 



Theorem 2.9. Under assumptions (M1)-(M5), suppose that condition (C) 
holds, the 2x2 matrix (a*,a 2 ) is nonsingular and (bi,b 2 ) 7^ (0,0). Then 
for some open interval I 3 there exist differentiable functions <pj : I — > R 
u»& 0j(O) = 0, j = 1 or 2, and : I -> M 2 u»& 0(0) = 0, 0"(O) ^ 
and 0(a) = (f)(—a) for a G I, such that a homoclinic orbit exists on X + for 
Hi — 01(^2) or [i 2 = 2 (/Wi); and an S-conjugate pair of homoclinic orbits 
exist for fi = 0(a): a pitchfork bifurcation of homoclinic orbits occurs (see 
Fig. 0). 
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a 

Figure 4: Pitchfork bifurcation of homoclinic orbits 

Proof. As in the proof of Theorem 12.71 we use Lemma 12.21 to obtain 
D M (id-n)F(^ ;0;0) 

= ^D^ - B x f(x h (t); 0)D„z - (id - U)DJ(x h (t); 0) = 0, 
D2(id-n)F(^ ;0;0) 

= 1d^ - DJ(x h (t); 0)D 2 a z - |(id - U)D 2 J(x h (t); 0)fo> 2 (t), ^(f)) = 0, 

where z = z(0, 0). Hence, we have (D^Zq)/-* — £ M an d L) 2 z — £ a - Noting 
that D a 2 = and using Lemma l2~2l we compute (ITTj) as 

Fj{ot,ii) 

/oo 
«,+;(*), -D M /(x h (t); 0)/i - cd^D,/^); 0)(//, p 2 (t)) 
•00 

- i« 2 D^(x h (t);0)(^ 2 (t)^ 2 (t)) - ia 3D '/(^W;0)(^(t),^ 2 (t),^(t)) 

- aD 2 /(x h (t); 0)(D^(t)/i + D^(t)« 2 , <^ 2 (t))>dt + G^o? + |/i| 4 ), 
so that 

/i) = ai/i + 6ia 2 + 0(^/a 6 + |/x| 3 ), 

F 2 (a,fi) = a 2 an + b 2 a 3 + 0( >/a 8 + |/i| 4 ). 

Applying Theorem 12.41 we see that a unique homoclinic orbits exists near 
(a, /x) satisfying 



a = 0, ai/i = or /j = — 
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Note that if there exists a homoclinic orbit which is not included in X + , then 
an S'-conjugate homoclinic orbit must exist. Thus, we repeat arguments given 
in the proofs of Theorems 12.51 and 12.71 to obtain the result. □ 



Remark 2.10. From Theorems \2. 7| and \2.9i we see that if condition (C) 
holds for n = 4, then a saddle-node or pitchfork bifurcation occurs under 
some nondegenerate condition, since uq = 2, as stated in Section 1. 

Remark 2.11. Suppose that ai/i + b\o? = 0. Then 

G{t; /i, a) = D M /(x h (t); 0)/i + \a 2 D 2 J{x\t); 0)(<p 2 {t), <p 2 {t)) 
satisfies UG(t; /i, a) = 0, so that £(t) = ^{t) + a 2 ^ a {t) is represented as 



m = m 



^*(r)G(r; /i, a) dr - a) 



where a) = a), . . . , a)) is given by 

/ (i>j(t),G(t; fM,a)) dt for j = 3, . . . ,n s ; 



Thus, we have 
a 2 H + b 2 a 2 



(tpj(t),G(t;n,a)) dt for j = n s + 3, . . . ,n; 



otherwise. 



(^ ns+2 (t),D M D x /(x h (t);0)(/i,^ 2 (t)) 

> 

+ la 2 D 3 x f(x h (t);0)(Mt),Mt),Mt)) 

+ D 2 j(x\ty,o)(at),Mt)))dt. 



(21) 



2.4 Hamiltonian case 



We finally consider the Hamiltonian case and set n = 2n, n £ N, and m = 1. 
We assume the following. 

(M6) There exists an analytic function H : M n xl->K such that 

f(x;iM) = J n (D x H(x;iM))*, 
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where J is the n x n symplectic matrix 



Jn 



id f 
-id s 



Assumption (M6) means that Eq. ([T]) becomes a Hamiltonian system 

x= J n (B x H(x;^Y (22) 

with a Hamiltonian function H(x; fi). Under this assumption we have n s = 
n n = n. See, e.g., [20] for more details on Hamiltonian systems. The n- 
dimensional stable and unstable manifolds, H /s (0) and H^ u (0), generically 
intersect along a one-dimensional curve, i.e., a homoclinic orbit, since they 
are included in an (n — l)-dimensional level set {x € lR n | H(x; /i) = H(0; fi)}. 

For the Hamiltonian system ( 122]) . the VE ()2]) and its adjoint equation ([6]) 
become 

t= J n DlH(x h (t);0)Z (23) 

and 

Z = D 2 x H(x\ty,0)U, (24) 

respectively, since T) 2 x H[x] 0) is a symmetric matrix and J* = —J n . Using 
solutions of the VE (|23|) . we can easily obtain solutions of the adjoint equation 
(El) as follows. 

Lemma 2.12. Let £ = ip(t) be a solution of ( )23|) . Then £ = J n <fi(t) is a 
solution of ( )24|) . 

Proof. Since £ = <p(t) satisfies (|23|) . we have 

= J n D^(x h (t);0)(^(t). 

We use the equality J 2 = — id n to modify the above equation as 

J n <p{t) = B 2 x H(x h (t); 0)J n (J n ip(t)) , 

which means that J n <p{t) is a solution of fl2l]). □ 
It is clear that these solutions are orthogonal 

(<p(t),J n <p{t)) = 0. 
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Hence, we take 

i/jn+j(t) = -J n <Pj(t), j = l,...,n , (25) 

and in particular 

*/j n+1 (t) = D x H(x h (t);Oy (26) 

since ip r (t) = x h {t) = J n B x H{x h {t); 0)*. 

We now turn to the problem of bifurcations in f )22|) . Introduce a new 
parameter v\ to modify f l2"2"j) as 

i = J n D x H(x;n)* + i/ 1 D x H(x;fj)*, (27) 

which depends on a two-dimensional parameter vector v = (fi, fi) G M 2 . 

Lemma 2.13. For 7^ 0, Eg. (j2"7|) /ias no homoclinic orbit passing through 
the set {x G W 1 \ D x H(x; /x) 7^ 0} ; especially in a neighborhood of x h (t) near 
fi = 0. 

Proof. Assume that v\ 7^ and let x(t) denote a homoclinic orbit in (|27|) . 
We easily compute 

^-H{x(t);fj,) = i^(D x H(x(ty,ii),T> x H(x(ty,ti)), 
at 

which is negative or positive for V\ < or > if D x H(x; /x) 7^ 0. This 
contradicts the fact that 

lim H(x(t);fi)= lim H(x(t);fi) 

since x(t) is a homoclinic orbit. Thus we obtain the result. □ 

By Lemma [2.131 no homoclinic orbit can be born from x h (t) for v\ 7^ 0. 
Now we apply Theorems O and EH to (ETJ). We use (J25D to have 

J (B x H(x h (t); 0), B x H(x h (t); 0))dt, 

and 

/"OO 

6 X =| / (D x H(x\ty0),Dl(J n D x H(x h (ty0))(Mt)iMt)))dt 



00 

00 



■00 

00 d 



00 dt 



J n D x H(x h (t); 0), D 2 ff (x h (t); 0)(^ 2 (t), p 2 (t))) dt 
(^ 2 (t),D 2 if(x h (t);0)^ 2 (t))dt = 
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since 



D^(x h (t);O)(0 2 (t),^ 2 (t)) 
= <0 2 (t),D^(x h (t);O)^(t)> 
= (ADlHix^O^^^lHix^O)^)) = 0. 



Using Lemma 12.131 and noting that the hypotheses of these theorems hold 
only if v\ = 0, we obtain the following theorems for (l22~l) with m — 1. 

Theorem 2.14. Under assumptions (M1)-(M3) and (M6), let n = 2 and 
suppose that a 2 ,6 2 ^ 0. Then a saddle-node bifurcation of homoclinic orbits 
occurs at fi = 0. Moreover, it is supercritical or subcritical, depending on 
whether a 2 6 2 < or > 0. 

Theorem 2.15. Under assumptions (M1)-(M6), letn^ = 2 and suppose that 
a 2 ,6 2 7^ 0. Then a pitchfork bifurcation of homoclinic orbits occurs at \i — 0. 
Moreover, it is supercritical or subcritical, depending on whether a 2 6 2 < or 



In these statements the constants a 2 ,fr 2 ,a 2 ,& 2 are obtained by ( fl9l) and 
(ED for ((22)) with m = 1 as follows: 



> 0. 




/ (J n¥ > 2 (t), D2/(x h (t); 0)(Mt), <P2(t))) dt, 



J — oo 




(28) 



+ Dlf(x h (t);Q)(e(t),ip 2 (t)))dt, 




+ D 2 J(x h (t);0)(C(t),Mt)))^ 



where we used (12!)|) with j = 2, and 




(29) 
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with G^t) = DJ(x h (t);0), G a (t) = iD2/(x h (t);0)(^ 2 (t),^ 2 (t)) and 

oo 

(i/)j(t), G^ait)) dt for j = 3, . . . , n s ; 

> 

(ipj{t), G^ a (t)) dt for j = n s + 3, . . . , n; 











x otherwise. 

See Remark 12.111 Note that in this case the condition a\V + b\0? = is 
equivalent to v\ = 0. 

3 Differential Galois theory 

Differential Galois theory deals with the problem of integrability by quadra- 
tures for differential equations. Here we briefly review a part of the differen- 
tial Galois theory which is often referred to as the Picard-Vessiot theory and 
gives a complete framework about the integrability by quadratures of linear 
differential equations with variable coefficients. 

3.1 Picard-Vessiot extensions 

Consider a system of abstract differential equations, 

y = Ay, A 6 gl(n,K), (30) 

where K is a differential field and gl(n, IK) denotes the ring of n x n matrices 
with entries in IK. We recall that a differential field is a field endowed with a 
derivation d, which is an additive endomorphism satisfying the Leibniz rule. 
By abuse of notation we write y instead of dy. The set Cr of elements of 
K for which d vanishes is a subfield of IK and called the field of constants of 
K. In our application of the theory in this paper, the differential field IK is 
the field of meromorphic functions on a Riemann surface V endowed with a 
meromorphic vector field, so that the field of constants becomes the field of 
complex numbers C. A differential field extension L D IK is a field extension 
such that L is also a differential field and the derivations on L and IK coincide 
on X. 

A differential field extension L D IK satisfying the following conditions is 
called a Picard-Vessiot extension for (]30p . 
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(PV1) There is a fundamental matrix $ of (I3"U|) with coefficients in L. 
(PV2) The field L is spanned by K and entries of the fundamental matrix $. 
(PV3) The fields of constants for L and K coincide. 

The system ( 1301) admits a Picard-Vessiot extension which is unique up 
to isomorphism. An algebraic construction of the Picard-Vessiot extension 
was given in a general situation by Kolchin (see, e.g., [16]). However, when 
K is the field of meromorphic functions, such a construction is unnecessary. 
Actually, we have a fundamental matrix in some field of convergent Laurent 
series, and get the Picard-Vessiot extension just by adding convergent Lau- 
rent series to the base field K. This is the point of view originally adopted 
by Picard and Vessiot (cf. [T3]). 

We now fix a Picard-Vessiot extension L D K and fundamental matrix 
$ with coefficients in L for fl30l) . Let o be a K- automorphism of L, which 
is a field automorphism of L that commutes with the derivation of L and 
leaves IK pointwise fixed. Obviously, cr($) is also a fundamental matrix of 
(130]) and consequently there is a matrix M a with constant entries such that 
cr($) = $M CT . This relation gives a faithful representation of the group of 
K-automorphisms of L on the general linear group as 

R: Aut K (L) -> GL(n, C L ), a ^ M a , 

where GL(n, Cl) is the group of n x n invertible matrices with entries in Cl- 
The image of R is a linear algebraic subgroup of GL(n, Cl), which is called 
the differential Galois group of (13"U|) and denoted by Gal(L/K). This repre- 
sentation is not unique and depends on the choice of the fundamental matrix 
$, but a different fundamental matrix only gives rise to a conjugated repre- 
sentation. Thus, the differential Galois group is unique up to conjugation as 
an algebraic subgroup of the general linear group. 

Definition 3.1. A differential field extension L D K is called 

(i) an integral extension if there exists a G L such that a G K and L = 
K(a) , where K(a) is the smallest extension o/K containing a; 

(ii) an exponential extension if there exists a G L such that a/ a G K and 
L = K(a); 
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(Hi) an algebraic extension if there exists a G L such that it is algebraic 
over K and L = K(a). 

Definition 3.2. A differential field extension LdIjs called a Liouvillian 
extension if it can be decomposed as a tower of extensions, 

L = K„ D . . . D Ki D K = X, 

swc/i f/iai eac/i extension Kj+i D Kj either integral, exponential or al- 
gebraic. It is called strictly Liouvillian if in the tower only integral and 
exponential extensions appear. 

Let G C GL(n, Cl) be an algebraic group. Then it contains a unique 
maximal connected algebraic subgroup G°, which is called the connected 
component of the identity or connected identity component. The connected 
identity component G° C G is a normal algebraic subgroup and the small- 
est subgroup of finite index, i.e., the quotient group G/G° is finite. By the 
Lie-Kolchin Theorem [131 EH] , a connected solvable linear algebraic group is 
triangularizable. Here a subgroup of GL(n, Cl) is said to be triangularizable 
if it is conjugated to a subgroup of the group of upper triangular matri- 
ces. The following theorem relates the solvability and triangularizability of 
the differential Galois group with the (strictly) Liouvillian Picard-Vessiot 
extension (see [131 EHj and [5] for the proofs of the first and second parts, 
respectively). 

Theorem 3.3. Let L D IK be a Picard-Vessiot extension of (I30p . 

(%) The connected identity component of the differential Galois group Gal(K/L) 
is solvable if and only ifhDWLisa Liouvillian extension. 

(ii) If the differential Galois group Gal(K/L) is triangularizable, then L D 
K is a strictly Liouvillian extension. 

3.2 Monodromy group and Fuchsian equations 

Let K be the field of meromorphic functions in a Riemann surface T and let 
to £ r be a nonsingular point in (1301) . We prolong the fundamental matrix 
analytically along any loop 7 based at to and containing no singular 
point, and obtain another fundamental matrix 7 * So there exists a 

constant nonsingular matrix M[ 7 ] such that 

7 * $(t) = $(t)M [7] . 
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The matrix M[ 7 ] depends on the homotopy class [7] of the loop 7 and is called 
the monodromy matrix of [7]. 

The set of singularities in (130]) is a discrete subset of T, which is denoted 
by ,y. Let ni(T \ ,y,t ) be the fundamental group of homotopy classes of 
loops based at to- We have a representation 

R: 7n(r\ y,t ) GL(n,C), [7] ^ M [7] . 

The image of i? is called the monodromy group of (1301) . As in the differential 
Galois group, the representation R depends on the choice of the fundamental 
matrix, but the monodromy group is defined as a group of matrices up to 
conjugation. In general, monodromy transformations define automorphisms 
of the corresponding Picard-Vessiot extension. 

A singular point of (130]) is called regular if the growth of solutions along 
any ray approaching the singular point is bounded by a meromorphic func- 
tion; otherwise it is called irregular. In particular, if A = B(t)/t with Bit) 
holomorphic at zero, then Eq. (13D]) has a regular singularity at t = 0. Eq. (1311 
is said to be Fuchsian if all singularities are regular. Any univalued solution 
of a Fuchsian equation is meromorphic. This gives us the following result 
along with the normality of the Picard-Vessiot extensions (see, e.g., Theo- 
rem 5.8 in [5U] for the proof). 

Theorem 3.4 (Schlessinger). Assume that Eq. (130]) is Fuchsian. Then the 
differential Galois group of (130]) is the Zariski closure of the monodromy 
group. 

Since the group of triangular matrices is algebraic, the Zariski closure of 
a triangularizable group is triangular izable. Noting this fact, we obtain the 
following result immediately from Theorem 13.41 

Corollary 3.5. Assume that Eq. (130]) is Fuchsian. Then the monodromy 
group is triangularizable if and only if the differential Galois group is trian- 
gularizable. 



4 Proof of Theorem 1.1 



We now prove our main result. Here we note that / and ^ can be extended 
to complex analytic ones in a neighborhood of M 4 in C 4 since they are real 
analytic. In this section we consider such complexifications including one of 
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Figure 5: Generators Uj, j — 1, 2, 3, 4. 



4.1 Normal and tangential variational equations 

We first use assumption (A2) to decompose the four-dimensional VE ([2]) into 
two-dimensional normal and tangent parts, so that we reduce our analysis to 
a two-dimensional system. 

Consider the generic variational equation at fi = 0, 

i = /O;0), i = T>J(x;0)£ t (31) 

which defines a flow on the tangent bundle TC 4 and which is linear on its 
fibers. Since ^# is invariant under the flow of ([T|) by assumption (A2), the 
tangent bundle and tangent bundle of C 4 restricted to TC 4 |^, 

are invariant under the flow of (13~T|) . The normal bundle is identified 

with the quotient TC 4 \^/T^ by definition, and its fiber at x G ^# is the 
linear space T X C 4 /T X ^ . Let us take a moving frame on i.e., a system 
of generators Uj G C 4 , j = 1, 2, 3, 4, for the tangent space T X C 4 with x G 
such that T x j& = span{ui,u 2 } and N x ^ = span{w 3 , w 4 } (see Fig. [5]). We 
introduce new coordinates (xi, X2, Vii V2) G C 4 by 

f = + X2U2 + ?7l^3 + ??2^4 

on TC 4 | ^. The invariance of T^t under the flow of fl3T|) ensures that the 
plane r\\ = 772 = is invariant in the restriction of ( 1311 to So the second 
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equation of (1311) is rewritten by a block form in the new coordinates as 



X2 



f A x {x) A c (x 




A v (x) 



fxx\ 

X2 

Vi 

W 



(32) 



if restricted to Here the 2x2 matrix functions A x (x), A v (x) and A c (x) 



are analytic and obtained via algebraic and differential manipulation from 
D x f(x; 0) and Uj, j = 1, 2, 3, 4. 

Let ej be a unit eigenvector of D x /(0; 0) corresponding to the eigenvalue 
Xj for j = 1, 2, 3, 4. From (1321) we obtain the following result. 



Lemma 4.1. The tangent space of ^ at the origin x = is spanned by e,- 



and ej + , i.e., Tq 



span{e 



3-1 °J+ 



}, where j„ = 1 or 2 an<i j_ 



3 or 4. 



Proof. Let x = in f[32|) . We easily see that two eigenvalues of the Jacobian 
matrix 

'A x (0) A c (0)' 
A v (0) / 

are those of A x (0), and the associated eigenvectors have the form £ = 
(x±; 0) G C 2 x C 2 , where x± are eigenvectors of A x (0). Since ^ contains a 
homoclinic orbit, A x (0) must have a pair of positive and negative eigenvalues. 
Thus we obtain the result. □ 

Henceforth we denote \± = A J± and write the other eigenvalues as fi±, 
where /i_ < < 

Let us consider (1321) for x = x h (t). Taking the normal components 77 = 
(V11V2), we have the normal variational equation (NVE) around x h (t), 



77 = A^ityr). 



(33) 



Moreover, we set \ = {X11X2) and V = to obtain the tangent variational 
equation (TVE) around x h (x), 



X = A x (x h (t)) X , 



(34) 



which governs the dynamics of fl3Tl) on for x = x h (t). When interested 
in a necessary condition for condition (C), we only have to deal with the 
two-dimensional NVE (121) instead of the four- dimensional VE (T51) as follows. 
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Lemma 4.2. If condition (C) holds, then 

(C) the NVE fl33|) has a non-vanishing bounded solution. 

Proof. Since A x (0) = \im t ^±oo A x (x h (t)) has two eigenvalue A_ < < A + 
(cf. the proof of Lemma l4TTj) . we see that the TVE (1341) has one independent 
bounded solution at most, as in Lemma 12. 1[ In addition, the solution £ = 
x h (t) of the VE © gives a bounded solution of the TVE (JMD- Hence, the 
TVE (I34p has only one independent bounded solution. 

Let us assume that condition (C) holds, i.e., the VE ([2]) has two inde- 
pendent bounded solutions. Since the TVE (]34j) has only one independent 
bounded solution, the NVE f[3"31 must have a non-vanishing bounded solu- 
tion. □ 



4.2 Analyses of the VE, NVE and TVE 

We next analyze the VE fl3J, NVE (PD and TVE flM} by using adequate 
systems of coordinates. We begin with the following lemma. 

Lemma 4.3. There exist analytical functions h± : U — > C 4 defined in a 
neighborhood U of x = in C stzc/i that h± (0) = and the complexification 
of the homoclinic orbit x h (t) is represented as 




h+ (e A +<) for Re(t) > 0; 
h- (e A - 4 ) /or Re(t) < 



Proof. Let T denote the complex separatrix of x = on the invariant mani- 
fold ^ in ([T]). It follows from assumption (A2) that x = is a double point 
on T. Hence, the intersection of T \ {0} with a small neighborhood of x = 
in C 4 consists of two pointed disks We add two different copies of the 
origin to these pointed disks, and get two disks E±. Using Lemma [4.11 and 
holomorphic changes of coordinates, we represent the restriction of (0Q) to E± 
as 

i ± = A ± z± + g±(z±)zl, (35) 

where z± represent coordinates in the disks E± and g±{z) are holomorphic 
functions. Since Re(A±) ^ 0, Eq. ( 1351) is analytically linearizable (see, e.g., 
Section 1.2 in Chapter 4 of [3]) and hence rewritten as 

C± = A±C± (36) 
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Figure 6: Local Riemann surface Fi oc and its covering {£±, E t }. 



under holomorphic changes of coordinates, by reducing the sizes of E± if 
necessary. The flow of (I36p is given by (±{t) = e x±t and the functions h± are 
the inverses of the transformation x h- > £± in U. □ 

Let r = {x = x h (t) 1 1 G M} U {0}. The curve T in the complex space 
C 4 consists of the homoclinic orbit x h (t) and saddle x — 0, which is a double 
point. We introduce two points + and 0~ corresponding to the origin for 
desingularizing the curve r . The points + and 0_ are represented in the 
temporal parametrization by t = +oo and t = — oo, respectively. Take a 
sufficiently narrow, simply connected neighborhood T, t of T \ (S + U S_) 
in T, such that it contains no singularity of ([2]) and intersects E± in two 
simply connected domains. We set ri oc = S_ US t US + , so that ri oc is simply 
connected and contains only two singularities of the VE (j2J) at 0±. See Fig. |6j 

Using Lemma [4.31 and the covering {E-j-,E t } of Fi oc , we introduce three 
charts on the Riemann surface Fi oc as follows: The charts in Y,± are given by 
z = e x±t and the chart in E t is given by t. Thus, we transform the VE (J2j) 
onto T\ oc , such that it is unchanged in E 4 and written as 



dz 



\±z± 



D x f{h±{z ± );0)£ 



(37) 



in E±. We easily see see that D x f(x±( 
following lemma. 



0) is analytic in z, and obtain the 



Lemma 4.4. The singularities of the VE (j2J) at 0± are regular. 



Thus, the VE 



is Fuchsian on Fi nr , and so are the NVE ([33]) and 



TVE ( 1341) . We also have the following result. 
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Lemma 4.5. By the coordinates z±, the NVE ( 13 3 p and TVE (T34"]) are, re- 
spectively, rewritten as 



— A, (*±)»7 — = — £ x (z±)x 



dz± z± v dz± z± 

in E± ; where the 2x2 matrix functions B^(z) and B^(z) are holomorphic 
and 

Proo/. As in ([ST]), we can rewrite the NVE ([33]) and TVE fl^U) as 

d.7? 1 d% 1 

-A v (h±(z±))r] and - — = A x (h±{z±))x, 



dz± X±z± dz± X±z± 

respectively, in £ ± . Noting that A v (0) and A x (0) have eigenvalues \± and 
fi±, we prove the result. □ 



4.3 Application of the differential Galois theory 

We apply the differential Galois theory of Section 3 to the NVE (133]) and 
VE ([2]) on the Riemann surface ri oc . To this end, we need an auxiliary result 
for two-dimensional systems of the form 

p- = lB(z)y, yeC 2 , (38) 
dz z 

where the 2x2 matrix function B(z) is holomorphic at z — 0. The system 
( 138]) has a regular singular point at zero and its fundamental matrix near 
z = is expressed as 

= Y(z)z E , (39) 

where the 2x2 matrix function Y(z) is meromorphic near z = and E 
is a constant matrix. Using this expression, we can easily compute the 
monodromy matrix corresponding to an infinitesimal loop around z = 
as M = exp(2iiiE) (see [I], page 8). 

Lemma 4.6. Assume that B(0) has two real eigenvalues p± such that p_ < 
< p + . Then the following statements hold: 
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(i) Eq. (138 p has a solution y{z) which is bounded along any ray approaching 
z = and unique up to constant factors. 

(ii) Any other independent solutions of (|38p are unbounded along any ray 
approaching z — 0. 

(Hi) The monodromy matrix M has an eigenvalue e 2mp+ , for which the as- 
sociated eigenvector is given by y{z), 

My(z)=e 2m ^y(z). 

Proof. We first consider the nonresonant case of p + — p_ ^ Z. Using Theo- 
rem 5 in Chapter 2 of [I], we have E = B(0) in fl39|) . Under a constant linear 
transformation, we assume that -B(O) is diagonal and 

*{z)=Y(z)(^ z ° p+ y (40) 

Letting $(z) = (2/1(2), 2/2(2)) and Y(z) = (vi(z), v 2 (z)) in (j4"U|) . we have 
2/1(2) = z p ~fi(z) and 2/2(2) = z p+ vi{z)- Since 1*1(2), 1*2(2) are holomorphic 
and z p+ (resp. z p ~) is bounded (resp. unbounded) along any ray approaching 
z = 0, so is the solution 2/2(2) (resp. 2/1(2))- Thus, we prove parts (i) and (ii). 
Moreover, part (iii) immediately follows by noting that M = exp(2niB(0)) 
is also diagonal in the present coordinates and y(z) = 2/2(2). 

We next consider the resonant case of p + — p_ G Z. Let po £ [0, 1) be a 
real number such that p± = po ± m-t for some positive integers m±. Using 
Theorem 6 in Chapter 2 of 0], we have a fundamental matrix of the form 

*{z) = Y{z)(^~ ^z E ' (41) 

under a constant linear transformation, where E' is a 2 x 2 matrix of the 
Jordan form having a double eigenvalue p. If £" is diagonal, then the expres- 
sion (l4"Tj) coincides with (l4"0"j) . so that we can prove parts (i), (ii) and (iii) as 
above. 

Suppose that E' is nondiagonal. This situation corresponds to so-called 
logarithmic singularity. Eq. (14 ip becomes 

*M = y MG4og; X)' < 42 > 
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and the same argument as in the nonresonant case yields parts (i) and (ii). 
Moreover, the monodromy matrix is computed as 



_ f27Ttp \ _ / e 2 ^° \ 
m - exp ^ 2m - \^ 2lxm 2,i P0 e 2,i P0 J ■ 

Letting Y(z) = (^i(z), w 2 (z)) as above, we see that 2/2(2) = e 2mpo V2(z) is an 
eigenvector of M for the eigenvalue e 2mp0 = e 2nlp+ . Thus we complete the 
proof. □ 

We turn to the NVE (1331) and VE (TSl) on the Riemann surface T 



loc- 

Theorem 4.7. Under condition (C), the monodromy group of the trans- 
formed NVE ( |33l) on Ti oc is triangularizable. 

Proof. Suppose that condition (C) holds, i.e., the NVE (133]) has a non- 
vanishing bounded solution. Using Lemmas 14.51 and 14.61 we see that the 
corresponding bounded solution of the NVE (133]) is a common eigenvector of 
the monodromy matrices around 0±. Since a group of 2 x 2 matrices with 
a common eigenvector is triangularizable, so is the monodromy group of the 
NVE fl3l on T loc . □ 



Proof of Theorem 11.11 

Suppose that condition (C) holds. Then it follows from Lemma 14.21 that 
condition (C) holds. Using Corollary 13.51 and Theorem 14 .7\ we see that 
the differential Galois group of the NVE (133]) on ri oc is triangularizable. 
Recall that the TVE ( 1341) is always integrable and solutions of the VE §2§ 
are obtained by substituting solutions of the NVE ( 13"3"]) and by solving the 
resulting equation by variation of constants. Hence, the differential Galois 
group of the VE ([2]) on ri oc is triangularizable. This completes the proof. □ 

Remark 4.8. Suppose that the full Riemann surface T in C 4 has genus zero, 
and let T be its desingularization. Then T is a Riemann sphere and hence 
the field of meromorphic functions on T is isomorphic to that of rational 
functions C(z). Assume that the pullback of the VE (T5]) to T has exactly 
three regular singularities {0±,p}. Then T\ oc \ {0±} is homotopic to T \ 
{0±,p}, so that both Riemann surfaces give rise to equivalent monodromy 
representations. Hence, we see via Theorems\l.l\ and\3.4\ that when regarded 



as a differential equation with coefficients in C(z), the VE (J2J) is integrable by 
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Liouvillian functions on C(z), which lie in a Liouvillian extension ofC(z), 
if condition (C) holds. This situation happens in the example given in the 
next section. 



5 Example 

To illustrate our theory, we consider 

xi = x 3 , x 3 = xx- (x\ + Pix\)x x - f3 3 x 2 , ^ 
±2 = x 4 , x A = sx 2 - {Pixl + P 2 x\)x 2 - /3 3 Xi - /3 4 X2, 

which represents steady states in coupled real Ginzburg-Landau equations of 
the form 



d t Ut = dlU x -C7i + {Ul + PiU\)U x + 3 U 2 , 
rd t U 2 = d 2 x U 2 - sU 2 + {(5^1 + P 2 Ut)U 2 + + p A U. 



2 

2 ■ 



U U U 2 ER, (44) 



with xi = Ui, x 2 = U 2 , x 3 = d x Ui and x 4 = d x U 2 , where s, r and fij, j = 
1, 2, 3, 4, are constants. See [2T] and references therein for general information 
on Ginzburg-Landau equations and [U [26] for examples of PDEs of the type 
( |44"|) . Eq. ( |4"3]) also represents a two-mode truncation of non-planar vibrations 
of a buckled beam when (3 3 = @4 = (see, e.g., [32]). Specifically, we assume 
that s > 1 since the case of s < 1 can be treated very similarly. 
Eq. ( I4"3"j) is Hamiltonian with a Hamiltonian function 

H = \{—x\ - SX\ + Pix\x\ + x\ + X4) + |(xf + /3 2 ^2) + /?3^1^2 + |/?4^2- 

The eigenvalues of the Jacobian matrix of the right-hand-side of (T43!) at x = 
are given by 

Ai = -y/s, X 2 = -1, A 3 = 1, A 4 = v^- 

When j3 3 = 0, the (xi,x 3 )- and the (x 2 , x 4 )-planes are invariant under the 
flow of fT43|) and there exist a pair of homoclinic orbits 

x\ {t) = [±V2 secht, 0, sech t tanh t, 0) (45) 

on the (xi, X3)-plane. Here we are only interested in the homoclinic orbits 
given by (j4"5"l) although another pair of homoclinic orbits exist on the (£2, X4)- 
plane. Thus, assumptions (Al) and (A2) (i.e., (Ml) and (M2)) hold for 
(3 3 = as well as assumption (M6). 
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Let p 3 = /3 4 = 0. Then Eq. 



,5' 



(\ 





is Z 2 -equivalent with 






-1 
1 




0\ 



-V 



Thus, assumption (M3) holds. We also have X + = {(xi, 0, £3, 0) G M 4 | x±, £3 G 
K}, X- = {(0, x 2 , 0, x 4 ) G 1R 4 I x 2 , x 4 G R}, and a£(t) G X+. Henceforth we 
take x h (t) = x\(t) to avoid complication in notation. 

For /3 3 = /5 4 = 0, the VE around the homoclinic orbits x h (t) is given by 

£i=£ 3 , e 3 = (l-6sech 2 t)£i, (46a) 
6 = e4 = (s-2/3 lS ech 2 t)e 2 , (46b) 
which is divided into uncoupled two-dimensional linear systems of the form 
Vi = V2, V2 = (vi - v 2 sech 2 t)?7i. (47) 

Eqs. fl46aD and fl46b|) . respectively, give the TVE and NVE for Let- 
ting 2 = sech 2 t in (1471) . we have a Fuchsian second-order differential equation 

d 2 ?7 



3z — 2 dr] V\ — v 2 z 
dz 2 + 2z{z-\)Tz + 4z 2 {z-l) V 



0. 



(48) 



which has regular singularities at z = 0, 1, 00 and whose solutions are ex- 
pressed by a Riemann P function [31] as 







1 





00 



,+ 1 



where Sq and represent the local exponents of ( 1481) at z 
respectively, and are given by 



(49) 



and 00, 



±1 



1 ± V4z/ 2 + 1 



Note that the local exponents of (f48l) at z — 1 are and ^. Thus, the 
VE (I46al b) is transformed into the Fuchsian system of the type (T48|) with 
three regular singularities, so that we only have to discuss the monodromy 
and differential Galois groups of (j48|) on the Riemann sphere C U {00} (see 
Remark 14. 8p . The following result was essentially proven in 
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20 



4 - 




s 



Figure 7: Condition (I5TI) for £ =0-4. The value of £ is labeled. Saddle-node 
bifurcations can occur on the red curves but not on the blue curves while 
pitchfork bifurcations can occur on all curves. 



Lemma 5.1. Consider a Fuchsian second-order differential equation which 
has three singularities z = Zj, j = 1, 2, 3, and a Riemann P function 




Then its monodromy and differential Galois groups are triangularizable if 
and only if at least one of pi + p 2 + p%, —p\ + p 2 + P3, Pi — P2 + P3 o,nd 
Pi + P2 — Pz is an °dd integer, where pj = pf — Pj, j = 1,2,3, denote the 
exponent differences. 



Using Lemma 15. 1[ we see that the monodromy and differential Galois 
groups for PH|) are triangularizable if and only if at least one of 



ISA - S 



-»o;±(*£-0±5= o ( 5 °) 



is an odd integer. Obviously, this condition is satisfied by (146 aft , and also by 
( I46bl) if and only if for some £ e Z 



V8/3i + 1 ± 2v^ ± 1 = 2{2£ + I) 
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ft = M±gtl)lzl , (eZ . (5i) 

See Fig. [7] for how the values of fli satisfying condition (1511) with £ =0-4 
change when s is varied. Noting that Eq. (CD) is Z2-equivalent for /?3 = and 
applying Theorem 11.11 we prove the following theorem. 

Theorem 5.2. [(%)] 

1. Choose /?3 as a control parameter and fix the other parameters. Then a 
saddle-node bifurcation occurs at /3 3 = in fj4*3|) on/?/ z/ condition floTl) 
/io/c/s. 

Choose 0i as a control parameter and fix the other parameters, espe- 
cially /?3 = /?4 = 0. T/ien a pitchfork bifurcation occurs in (j4*3"]) on/?/ i/ 
condition f[5"Tj) holds. 

Now we compute the second independent bounded solution of the VE 
(I46al b) for /3 3 = /3 4 = 0, based on the above result. Consider ( 146b[) and 
suppose that condition ( 15T|) holds. Then we have 

4 = ±±v^, 4 = i (1 ± |2v^ + 2£ + 1|) 

in Eq. §DJ- We set 

( = z s or], z s o (z - l) 1/2 n, z s or] and z s ° (z - l) 1/2 n, 

so that the Riemann P function fT49]) becomes 



z s ° P < s 



z s ° P { s 
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and 



1 oo 1 f 1 oo 



respectively. Hence, Eq. ( |48|) is transformed to the hypergeometric equation 

d 2 C dC 
z(l - + (c - (a + 6 + l)z)-g - a&C = 0, 

where c = ±(sq — Sq ) + 1 = 1 and 

(a, 6) = + 1), -Vi - \i) , - 1), V~s + |* + 1) , 

(-i£,v^+ + and (|£ + l,-v^-i(£-l)), 

respectively. Using a well-known result on the hypergeometric equation (e.g., 
[3~T]). we obtain a bounded solution of (I46b|) as 

f 2 (i) = ^/ 2 (1 - 2 )V2 F (_A; + l, ^ + A; + J, 1; (52) 

for £ = 2 A; — 1 and 

^)=2^ /2 F(-Hl,^+fc-i,M) (53) 

for I = 2(k — 1), where k £ N, z = sech 2 t and F(a, b, c; z) is the hypergeo- 
metric function, 

^ ' £S j!c(c+1)...(c + j-1) 

which becomes a finite series when a is a nonpositive integer. Note that 
Eq. (|48p also allows a solution of finite series as 



&(*) = z~^ 2 F(-k + 1, -a/s + fc — |, 1; z) 
for £ = -2fc + 1 and 

for £ = —2k with k £ N but they are unbounded since z _1 = cosh 2 t — > oo 
as t — > 0. Thus, for /3 3 = /3 4 = 0, if fii satisfies flBTl) with £ > 0, then 
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condition (C) holds and the second independent bounded solution of the VE 
(jjgaj b) is given in the form <p 2 (t) = (0, f 2 (i), 0, &(*))* b Y §3) or 

We next carry out the Melnikov analysis for (I43p . Obviously, assump- 
tions (M1)-(M3) and (M6) hold, and assumptions (M4) and (A5) for fa = 
fa = 0. Using (EHD, we have 

/OO POO 
&(t)4(t)dt, b 2 = -fa / g(t)dt (54) 
-oo J —oo 

if we take /i = fa as a control parameter, where £2^) is given by ( 152]) or 
f )53|) . Since x\(t) is an even function of t and £ 2 (t) is an even or odd function 
depending on whether I (> 0) in ( l5"Tj) is even or odd, we easily see that 
a 2 , &2 7^ only if I is even and /?4 7^ 0. See Appendix A for computations of 
a 2 , 6 2 when £ = 0,2, 4. 

On the other hand, assume that fa = fa = and let /j, — fa be a control 
parameter. We take 

(fl(t) =(secht tanht, 0, 2 sech 3 t — secht, 0), 

(Ps(t) = (|t sech t tanh t + | sinh t tanh t — secht, 0, 

3t sech 3 t + 3 sech t tanh t — |t sech t + | sinh t) 
^1 (t) = (3i sech 3 t + 3 sech t tanh t — |t sech t + | sinht, 0, 

— |t sech i tanh i — | sinh i tanh t + sech t, 0) , 
ijj*(t) = (-2sech 3 t + sech t, 0, sech t tanh t, 0) 

(see [9]) so that 

ff(t) = M*) /Vi3(r)x^(r)£ 2 2 (r)dr + ^ 3 i(t) /^(^(r^dr, 
Jo Jo 

where and ^ (t) are the j-th components of y?i(t) and ^(Oj respec- 

tively, while = since D /J /(x h (t); 0) = 0. Hence, 

/oo 
e 2 2 (t) [*?(t)] 2 dt, 

00 ('551 

/•oo /-oo yuu j 

b 2 = -2 fa / x\(t)g(t)£(t)dt - fa / tf(t)dt. 



Note that £f(t) is an even function of t since (pu(t), ^33 (^) are odd and 
V3i(^)) ^13 (^) are even. We easily see that a 2 < and fr 2 7^ for almost all 
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values of 2 although the integral including in (|5"5"|) is difficult to be 

estimated analytically. 

Applying Theorems 12.141 and 12.151 we obtain the following theorem. 

Theorem 5.3. [(%)] 

1. Choose 3 as a control parameter and fix the other parameters. Then 
a saddle-node bifurcation occurs at 3 = in (143)) with 04 ^ if 
condition f[5"Tj) holds for a nonnegative and even integer I and a 2 ,b 2 ^ 
0. Moreover, it is supercritical or subcritical, depending on whether 
a 2 b 2 < or > 0. 

2. Choose 0i as a control parameter and fix the other parameters, espe- 
cially /?3 = 04 = 0. Then a pitchfork bifurcation occurs in (J32D with 
02 7^ if condition (l5Tj) holds for a nonnegative integer £ and b 2 ^ 0. 
Moreover, it is supercritical or subcritical, depending on whether b 2 > 
or < 0. 

In particular, if £ = or if i = 2, 4 and s is sufficiently large (s > 
0.16049 or s > 5.17784 x 10~ 7 ), then the quantity a2&2 has the same sign 
as 04 as shown in Appendix A, and the saddle-node bifurcations detected 
by Theorem I5.3( i) are subcritical (resp. supercritical) for 04 > (resp. for 

04 < 0). 

We finally give numerical results for fj4^|) . We used the computation tool 
AUT097 with HomCont [7] and performed continuations of homoclinic orbits 
with two parameters in a general setting for 

Xl = X-3, X 3 =Xi- (x\ + 0lX 2 2 )xi - 02,X 2 - ViX 3 , 

X 2 = X 4 , ±4 = SX 2 - {0ix\ + 2 x\)x 2 - 3 X! - 04x1 - U1X4 

instead of a system of the type ff27]) . Note that a statement similar to that 
of Lemma 12.131 still hold in (156]) and the homoclinic orbit persists only if 
V\ = 0. The homoclinic orbit fT45]) was taken as the starting solution in these 
continuations. Henceforth we fix the parameter s = 2. 

Figure [8] shows branches of homoclinic orbits when 0% is varied as a con- 
trol parameter and 0\ satisfies condition (!5Tl) for £ = 0,2 and 4. In plate (c) 
of Fig. [HI the maximum and minimum of x 2 (t) are plotted since x 2 (t) has 
no maximum and minimum at t — 0. From Fig. [H] we see that saddle- 
node bifurcations occur at 03 = 0, as predicted in Theorem 15.31 Moreover, 
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-0.2 

Figure 8: Bifurcation diagram with 3 a control parameter for s = 2: (a) 
/?! = 1.7071068, /3 2 = 1 and /3 4 = 2; (b) (3 1 = 7.5355339, /3 2 = 1 and /3 4 = 2; 
(c) ft = 17.36396103, (3 2 = 10 and /3 4 = 20. Note that condition (|5T|) 
approximately holds in plates (a), (b) and (c) for £ = 0, 2 and 4, respectively. 
In plates (c) two extrema varying continuously with (3\ are plotted in red and 
blue. 



these bifurcations are subcritical, as predicted by Theorem I5.3( i) with the 
computations of a 2 and b 2 in Appendix A. Note that x 2 (t) = at the bifur- 
cation point. We also remark that no saddle- node bifurcation was observed 
at j3 3 = when /3 4 satisfies condition (I5T|) with I = 1,3, and that secondary 
saddle-node bifurcations occur very near /3 3 = and the branch shape be- 
comes very different when I is higher and (3 2 , are smaller (this is the reason 
why large values of (3 2 , fi± are taken for in plate (c) of Fig. [8]). Profiles of ho- 
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Figure 9: Profiles of homoclinic orbits on the branches in Fig. [HJ (al,2) 
/3 3 = 1; (bl,2) /3 3 = 0.5; and (cl,2), fa = 0.3. Other parameter values in 
plates (al,2), (bl,2) and (cl,2) are, respectively, the same as in plates (a), 
(b) and (c) of Fig. EJ 



moclinic orbits on the branches in Fig. [HJ are given in Fig. [9j Note that these 
homoclinic orbits are symmetric. Asymmetric homoclinic orbits were also 
born from homoclinic orbits on these branches at pitchfork bifurcations as 
well as at the secondary saddle-node bifurcations although no branches of 
asymmetric orbits are drawn in Fig. [HJ 

Figure [10] shows branches of homoclinic orbits when /?i is varied as a 
control parameter for /3 2 = 1 and P3 = Pa = 0. Note that there exist a 
branch of X2{= £4) = for all values of 0±, and a branch which is symmetric 
about max X2 = to each one in Fig. [TQ1 From Fig. [TO] we see that pitchfork 
bifurcations occur when 0i satisfies condition ( 151]) for i =0-4, as predicted in 
Theorem 15.31 The first three bifurcations are supercritical, and the rest two 
ones are subcritical. Profiles of homoclinic orbits on the branches in Fig. [10] 
are given in Fig. [11] Note that homoclinic orbits on all branches in Fig. [10] 
are symmetric, like the profiles plotted in Fig. [TT] 
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Figure 10: Bifurcation diagram with f3\ a control parameter for s = 2, ft = 1 
and ft = ft = 0. 
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Figure 11: Profiles of homoclinic orbits on the branches of £ =0-4 in Fig. HO] 
for ft = 15: (a) £ = 0; (b) £ = 1; (c) £ = 2; (d) £ = 3 and (e) £ = 4 (the 
lower branch). Other parameter values are the same as in Fig. [TUl 
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A Computations of 62 given by ( ED for £ = 

0,2,4 

We first note that 

sech a t dt = 2 1 , 57 

r(a) 

where T(z) is the Gamma function, 

POO 

F(z) = / t^e-'dt. 
Jo 

Using the formula (1571) . we can compute 62 as follows. 

A.l Case of I = 

We set k = 1 in ( 1531) to have 

&(t) = sech^t, 

so that 



a -2 



,00 9v /i + 1 /2r 2 ('l A /i+ I) 

-V2 sech^tdt = - 2 r( ^* +2 \ 

r°° 2 3 ^ _1 r 2 f- /s) 

b 2 = -fa / sech 3 ^ dt = — W ' - fa. 

J-00 r(3Vs) 

Hence, the quantity 02^2 has the same sign as fa since T(z) > for z > 0. 

A. 2 Case of £ = 2 

We set k = 2 in (E3J) to have 

= sech^t sech 2 t) , 
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so that 



a 2 = - f (sech^+H - (Vs + |) sech^ +3 t) dt 

2^+1/2^(1^ + i j 2^+5/2^(1^ + 1; 



+ (>/?+§) 



r(^ + i) vv 27 r(^ + 3) 
_ 2^-v 2 (2 g + 3y/i - i)r 2 (^ + 1) 
(Vi + 2)r(v^+i) 

6 2 = - /3 4 / (sech 3 ^ -3(v^+|) sech 3v/i+2 t 

+ 3 (Vi + D'sech 3 ^ - |) 3 sech 3 ^+ 6 t)dt 

/ gv^jyg 3 , 2 3 ^r 2 (| v ^+i) 
- _/ H r(3V5) ~ 3 ^ s + -^ r(3Vi + 2) 

,o//-, 3^2^51^+ 2) 3 3 2^1^+3) \ 

+ 3(Vs + 2j r(3^ + 4) Ws + * } r(3^ + 6) J 

_2 3 v /i " 4 (72s 3 + 252s 5 / 2 + 262s 2 + 93s 3 / 2 + 72s + 32^ - 40)r 2 (|v^) 
= (3^ + l)(v^ + l)(3v^ + 5)r(3v/i) ^ 

where we used the relation F(z + 1) = zT(z). We see that the quantity a^hi 
has the same sign as /3 4 for s > 0.16049. 

A. 3 Case of I = 4 

We set k = 3 in fl53l) to have 



&(t) = sech^t (1 - (2>/s + 5)sech 2 t + (Vs + |) (Vs + I) sech 4 t) 
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so that 
a 2 



/oo 
(sech^ 1 * - (2v^ + 5)sech^ +3 t 
■oo 

+ {yfs + \) (v^+ I) sech^ +5 t)dt 

2^+^2(1^+1) 9^+5/2^(1^+1) 



(2^ + 5)- 



r(Vi+i) v v ; r(v^ + 3) 

2 v^+9/2 r 2/l /-, 

(v^+|) W-8 + D- 



r(v^ + 5) 

_2^-3/2( 4s 2 + 48s 3/2 + ig9s + 320^ + 153)r 2 (iv^ + \) 
= ( v ^ + 2)( v ^ + 4)r( v ^+l) 

/oo 
sech 3v/i t(l - (2v^ + 5)sech 2 t 
-oo 

+ (Vs + f) (v^+I) sech 4 t) 3 dt 

/ 2V^r 2 (jyi) 2V^r 2 (§ ^ + 1) 

= _/ H r(3V5) 3 ^ + 5 )- 

+ (v^+|) (5v^+f) 



T(3 V / i + 2) 

2 3 ^+ 3 r 2 (|^i + 2) 



2(v^+|) 2 (l0v^ + 3l) 



r(3v^ + 4) 
2 3 ^+ 5 r 2 (|v^ + 3) 

r(3^i + 6) 

2 3 ^+ 7 r 2 (|Vi + 4) 



+3( v5 + r(vi + i)(v5 + f) r(3yi+8) 

2V^( v /s")r 2 (§ v /i) 



(3 v ^+l)( v ^+l)(3 v ^ + 5)(3 v ^+7)( v ^ + 3)(3 v ^+ll)r(3 v ^) 
where 

g(s) =5184s 12 + 176256s 11 + 2519568s 10 + 20488032s 9 + 106620652s 8 
+ 375344312s 7 + 915087795s 6 + 1546383098s 5 + 1772860056s 4 
+ 1308687720s 3 + 556461984s 2 + 102326688s - 73920. 



i- 
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We also see that the quantity 02^2 has the same sign as (3^ for s > 5.17784 x 
1(T 7 . 
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